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Luminous superclusters: remnants from inflation 
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Abstract. We derive the luminosity and multiplicity functions of superclusters compiled for the 2dF Galaxy Redshift Survey, 
the Sloan Digital Sky Survey (Data Release 4), and for three samples of simulated superclusters. We find for all supercluster 
samples Density Field (DF) clusters, which represent high-density peaks of the class of Abell clusters, and use median lumi- 
nosities/masses of richness class 1 DF-clusters to calculate relative luminosity/mass functions. We show that the fraction of 
very luminous (massive) superclusters in real samples is more than tenfolds greater than in simulated samples. Superclusters 
are generated by large-scale density perturbations which evolve very slowly. The absence of very luminous superclusters in 
simulations can be explained either by non-proper treatment of large-scale perturbations, or by some yet unknown processes in 
the very early Universe. 
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1. Introduction 

Superclusters are the largest density enhancements in the 
Universe of common origin. Superclusters evolve slowly and 
contain information from the very early Universe. The investi- 
gation of large systems of galaxies was pioneered by the study 
of the Local Supercluster by de Vaucouleurs (I1953I I. Another 
approach to define superclusters was initiated by Abell IT958I 
[1961), who considered them as "clusters of clusters". Until 
recently, superclusters have been found mostly on the basis 
of catalogues of rich clusters of galaxies by Abell ('1958') and 
Abell et al. (1989). Abell supercluster catalogues have been 
complied by Zucca et al. iT9^, Ein asto et al. (IT994I [T997bl 
UMT S and Kalinkov & Kuneva fl99?). 

Actually superclusters consist of galaxy systems of dif- 
ferent richness: single galaxies, galaxy groups and clusters, 
aligned to chains (filaments). This has akeady been realized be- 
fore by Joeveer, Einasto & Tago ( 1978) . Gregory & Thompson 
J1978> . Zeldovich, Einasto & Shandarin J1982> . and has been 
confirmed by recent studies of superclusters using new deep 
galaxy surveys, such as the Las Campanas Galaxy Redshift 
Survey, the 2 degree Field Galaxy Redshift Survey (2dFGRS, 
Colless et al. "200? '200T) and the Sloan Digital Sky Survey 
Data Release 4 (SDSS DR4, Adelman-McCarthy et al. 2006). 
New galaxy redshift surveys are almost complete in a fixed ap- 
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parent magnitude interval. This allows to estimate total lumi- 
nosities of superclusters using weights, inversely proportional 
to the number of galaxies in the observational window of ap- 
parent magnitudes. This possibility has been used in recent su- 
percluster studies by Basilakos J2003> . Basilakos et al. J2001> . 
Erdogdu et al. (2004), Porter and Raychaudhury ("20051). and 
Einasto et al. (2003a 2003b 2005 2006a hereafter Paper I). 

In Paper I we compiled a catalog of superclusters using 
the 2dFGRS. A similar catalogue on the basis of SDSS Data 
Release 4 (DR4) is in preparation (Einasto et al. 2006ci here- 
after E06c). The properties of 2dFGRS superclusters were an- 
alyzed by Einasto et al. (2006b Paper II). The parameters 
of 2dFGRS superclusters were compared with properties of 
model superclusters based on the Millennium Simulation of the 
evolution of the Universe by Springel et al. (2005 ). This com- 
parison of real superclusters with simulated ones shows that 
geometric properties of simulated superclusters agree very well 
with similar properties of real superclusters. This demonstrates 
that the ability to simulate processes which lead to the forma- 
tion of superclusters has reached an advanced stage. However, 
one property of model superclusters is in conflict with reality: 
real samples have many more very luminous superclusters than 
model samples. The presence of very massive superclusters in 
our neighborhood is well known, examples are the Shapley and 
Horologium-Reticulum Superclusters (see Fleenor et al. 1 20051 
Proust et al. l2006l Nichol et al. .2006. and Ragone et al. .2006l and 
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Fig. 1. The cumulative distribution of peak densities of DF-clusters. Left panel shows densities uncorrected for the relative bias, 
right panel densities corrected for differences in the relative bias. 



references therein). However, the number of such extremely 
massive superclusters was too small to make definite conclu- 
sions on the phenomenon. 

The goal of this Letter is to determine the luminosity and 
multiplicity functions of 2dFGRS and SDSS DR4 superclus- 
ters and to compare these functions with similar functions of 
simulated superclusters. These observational samples of super- 
clusters are the largest available today. To verify the robust- 
ness of the results obtained with the Millennium simulation 
data we also use superclusters derived from a cosmological 
simulation of the same volume but with a lower mass reso- 
lution. This large collection of real and model superclusters 
allows us to make definite conclusions on the statistics of lu- 
minous superclusters. Supercluster catalogues of the 2dFGRS, 
as well as fits files of luminosity density fields of 2dFGRS and 
Millennium Simulation are available electronically at the web- 
site http : //www . aai . ee/~inaret/2df scl . html. 



2. Data 

In this analysis we used galaxy and group samples of the 2dF- 
GRS (Tago et al. 2006a Paper I) and SDSS DR4 (Tago et al. 
I2006bl E06c), and three simulated samples. To compile super- 
cluster catalogues we used the density field method which al- 
lows to correct the field to take into account the presence of 
faint galaxies outside the observational window, for details see 
Paper I. The richness of a supercluster can be characterised by 
its total luminosity, and by the number of rich galaxy clusters 
in it, i.e. the multiplicity function of the supercluster. In Papers 
I and II we derived both these characteristics. The sample of 
Abell clusters covers only our close neighbourhood, in more 
distant superclusters they are found only in exceptional cases. 
To have a richness parameter we used instead of Abell clus- 
ters high-density peaks of the density field, which we call DF- 
clusters. The spatial density of DF-clusters is approximately 
close to the spatial density of Abell clusters, thus their utiliza- 
tion as richness indicator yields results comparable with previ- 
ous supercluster studies based on Abell clusters. 



For comparison to the observational data sets we used the 
galaxy catalogues based upon the Millennium Simulation, for 
details see Springel et al. ( I2005> . Gao et al. ( I2005t and Croton 
et al. (2006 1. In addition, we also performed a cosmological 
dark matter simulation of a computational volume of side- 
length 500 Mpc using 256^ particles (model M500). We 
chose the concordance cosmology with the parameters Q.„ - 
0.27, Qa = 0.73, erg - 0.84. The simulation was carried 
out with the open source Multi Level Adaptive Particle Mesh 
(MLAPM) code by Knebe et al. 02001 1 and DM-halos have 
been found by the conventional Friends-of-Friends (FoF) pro- 
cedure. We denote simulated samples as Mill.AS and M500. 
The density fields of these models were calculated using all 
simulation galaxies/particles. In addition, we used the model 
Mill.FS, where galaxies were chosen using similar selection 
criteria as in the 2dFGRS sample (see Paper I). 

To get comparable results for the luminosity and multiplic- 
ity functions of dififerent samples all data and data reduction 
procedures must be as similar as possible. First of all, density 
fields must have identical threshold bias levels. We use in our 
study relative densities expressed in units of the mean density 
of the particular sample. To check the relative bias levels we 
found for all samples the threshold density, which yields non- 
percolating superclusters of maximal diameter approximately 
100 - 120 ' Mpc, and selected DF-clusters using a threshold 
density about 10 % higher than used in the selection of super- 
clusters. 

The comparison of peak densities of DF-clusters in differ- 
ing samples used in the present study shows considerable dis- 
parities amongst the samples, see Fig. [2 These variations are 
due to differences in mean densities used in the calculation of 
relative densities, i.e. differences in the threshold bias factor 
(for a detailed discussion of this phenomenon see Einasto et 
al. 11999b . Fig.[2shows that peak density distributions of DF- 
clusters of the SDSS and Millennium simulation samples are 
practically identical, thus we have used these samples as a stan- 
dard. To bring the threshold biases of the other samples to a 
comparative level we divided the density fields of the 2dFGRS 
samples by 0.739, and the model M500 density field by 2.025 
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Sample 


V 


Ngal 






Do 
yv, 


= 5.0 


Id 




N,, 




A'l 


= 6.0 

Lo 


lid 


n.id 


2dFc 


30.3 


184395 


2555 


567 


279 


6.86e+ll 


84 


19 


1664 


544 


265 


8.20e+ll 


55 


18 


SDSS 


43.9 


368002 


3621 


1012 


517 


7.09e+ll 


82 


23 


2364 


911 


483 


8.11e+ll 


54 


21 


Mill.AS 


125.0 


8964936 


4914 


2259 


1292 


1.64e+12 


39 


18 


2878 


1733 


1025 


2.10e+12 


23 


14 


Mill.FS 


125.0 


2094187 


3020 


1299 


762 


3.08e+12 


24 


10 


1687 


1068 


752 


3.03e+12 


14 


8 


M500c 


125.0 


9785827 


3032 


1860 


1281 


4.75e+14 


24 


15 


1880 


1440 


1037 


5.81e+14 


15 


12 



(these values were found by comparing peak density distribu- 
tions shown in the left panel of Fig.[0. The distribution of peak 
densities of the corrected density fields is shown in the right 
panel of Fig.[l] We see that there are practically no systematic 
differences between various samples. These corrected density 
fields were used to compile the supercluster catalogues used 
throughout this study. 

To obtain comparable results we used for all samples iden- 
tical procedures in the preparation of the data. In all cases 
superclusters were found using a luminosity (or mass) den- 
sity field smoothed with an Epanechnikov kernel of radius 
8 /z"' Mpc. For observational samples densities were calculated 
using weights of galaxies which take into account galaxies and 
galaxy groups too faint to fall into the observational window 
of absolute magnitudes at the distance of the galaxy. This is 
the conventional approach for obtaining the luminosity density 
field (Basilakos et al. 2001 Paper I). The density field was 
found for a cell size of 1 /i ' Mpc, which allows to investi- 
gate in detail the internal structure of superclusters (see Paper 
II). Superclusters were defined as connected non-percolating 
systems with densities above a certain threshold density. This 
threshold density is similar to the linking length used in the 
Friends-of-Friends (FoF) method to find systems of galaxies 
(or particles in simulations). 

Supercluster catalogues have been compiled for two val- 
ues of the threshold density, i.e. Dq - 5.0 and 6.0 (in units 
of the corrected mean density). In both cases a lower limit of 
the supercluster volume of 100 (/;"' Mpc)^ was applied. The 
main characteristics of these catalogues have been summarized 
in Tabled the name of the sample (c indicates the use of cor- 
rected density fields as described above), V is the volume (in 
million cubic Mpc), Ngai is the number of galaxies (or par- 
ticles in case of the model M500c) used in the determination of 
the density field and the luminosity/mass of superclusters. Next 
we give for both threshold densities the number of DF-clusters 
in the sample, A^^.; (to derive DF-clusters we used threshold den- 
sities 5.5 and 6.5 for the supercluster samples 5.0 and 6.0, re- 
spectively). Nsd and A^i are the total number of superclusters, 
and the number of superclusters of multiplicity 1 (i.e. contain- 
ing just one DF-cluster), respectively. Lq is the mean luminos- 
ity (or mass for M500c) of superclusters of multiplicity 1, ex- 
pressed in Solar units. These luminosities/masses were used 
in the calibration of relative luminosity functions, see the next 
section for details. Finally, rid and Usd are the mean spatial den- 
sities of DF-clusters and superclusters: the number of objects 



per million cubic /; ' Mpc. The threshold density 6.0 for the 
corrected density field of the 2dFGRS sample corresponds to 
a threshold density of 4.4 in the uncorrected density field, very 
close to the value 4.6 used in Paper I to find 2dFGRS superclus- 
ters. Thus the number of superclusters and their parameters is 
very close to the respective data used in Paper I. 



3. The luminosity and multiplicity functions of 
superclusters 

We shall use in this paper two independent parameters to char- 
acterise quantitatively the richness of a superclusters: the mul- 
tiplicity and the total luminosity (or mass). We define the multi- 
plicity of a supercluster by the number of DF-clusters in it. DF- 
clusters are high-density peaks of the density field, smoothed 
on a scale of 8 /i ' Mpc. As seen from Table ^ for thresh- 
old density Do - 6.0 the spatial density of DF-clusters in 
our samples is about twice the spatial density of Abell clus- 
ters, 25 per million cubic Mpc (Einasto et al. I1997bt : for 
a threshold density of 5.0 the density is somewhat higher. The 
other integral parameter of a supercluster is its total luminos- 
ity or mass, determined by summing luminosities/masses of 
all galaxies and groups of galaxies (DM-particles) inside the 
threshold iso-density contour which was used in the definition 
of superclusters. The relationship between the multiplicity and 
total luminosity (or mass for M500) is presented in Fig.|2l We 
see that luminosities have a rather large spread for superclus- 
ters of given multiplicity, the lower the multiplicity the larger 
the spread. The other feature, seen in Fig.|2] is the shift in the 
mean luminosities of superclusters for different samples. Partly 
this may be caused by the use of different color systems in var- 
ious samples. 

We are interested in the relative fraction of rich and very 
rich superclusters in respect to the number of poor superclus- 
ters. To avoid complications due to the use of different color 
systems and masses in case of the model M500, we define rel- 
ative luminosities (masses) as the luminosity (mass) in terms 
of the mean luminosity (mass) of poor superclusters, i.e. super- 
clusters that only contain one DF-cluster and hence are classi- 
fied as richness class 1 . The distribution of luminosities is ap- 
proximately symmetrical on a logarithmic scale (see Fig. |5J. 
We therefore used the logarithm of the luminosity to derive 
the mean value. This value, Lo, is also listed for all samples 
in TableQ] in units of the Solar luminosity (mass). 
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Fig. 2. The relationship between the multiplicity and the total luminosity (mass) of real and simulated superclusters, calculated 
for threshold density 6.0. The samples are: 2dFGRS (top left), SDSS DR4 (top right). Millennium Simulation Mill.AS (bottom 
left), the DM simulation M500 (bottom right). 



In Fig.|3lwe now show the relative luminosity (mass) func- 
tions (left panels) alongside the multiplicity functions (right 
panels) for the observational and model samples. The spatial 
density of superclusters is expressed in terms of the total num- 
ber of superclusters in the respective sample to avoid small dif- 
ferences due to the mean number density of superclusters in 
different samples. In upper panels these functions are shown 
for all observational and model samples using threshold density 
5.0, in lower panels only for the combined observational sam- 
ples Obsl and Obs2 (see below), and model samples Mill.AS 
and Mill.FS, using threshold density 6.0. The comparison of 
data obtained with threshold densities 5.0 and 6.0 shows that in 
the first case some superclusters of the real sample are actually 
percolating (the maximal diameter of the largest superclusters 
is almost 200 Mpc). For this reason we have used in the fi- 
nal analysis supercluster samples found with threshold density 
6.0. 

Upper panels of Fig. |3l show that there are small differ- 
ences between the observational samples 2dFGRS and SDSS: 
the 2dFGRS sample contains relatively more rich and very rich 
superclusters. These differences may be due to the differing 
depth of these samples: the limiting magnitude of the 2dF- 
GRS sample is about 19.35, whereas that of the SDSS sam- 



ple is 17.7. Clusters form superclusters via intermediate den- 
sity galaxy bridges. The 2dFGRS sample contains more faint 
galaxy bridges between high-density regions which facilitates 
the formation of more luminous superclusters. A similar differ- 
ence is observed between the two model samples. The Mill.AS 
sample has a very high spatial resolution and contains numer- 
ous faint galaxies which also acts in favour of joining nearby 
high-density regions via galaxy bridges. 

To check this explanation of the difference between sam- 
ples of various depth we used the simulated 2dFGRS sample 
Mill.FS, containing only galaxies brighter than 19.35, and cal- 
culated for this sample luminosity and multiplicity functions as 
for other superclusters samples, using threshold densities 5.0 
and 6.0. These versions are also shown in Fig.|3] We see that 
in this case the fraction of luminous superclusters in lower than 
for the full sample Mill.AS. This test shows that faint galaxy 
bridges between clusters play indeed an important role in the 
formation of superclusters. To reduce the SDSS DR4 super- 
cluster sample to the same level of bridge strength as the 2dF- 
GRS sample, a lower threshold density must be used. Trial 
calculations showed that a combined observational sample of 
superclusters can be formed using the 2dFGRS sample with 
threshold density 6.0, and the SDSS DR4 sample with a thresh- 
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Fig. 3. Upper panels show the relative integrated luminosity/mass functions (left), and the multiplicity functions (right) of all real 
and simulated superclusters, found using threshold density 5.0. Luminosities/masses are expressed in units of mean luminosi- 
ties/masses of superclusters of multiplicity 1, Lq, given in Tabled Different lines mark supercluster samples: the 2dFGRS, the 
Sloan Digital Sky Survey DR4, the Millennium simulation full sample Mill.AS and selected sample Mill.FS, the M500 sample 
and the sample of Abell superclusters. Lower panels show the comparison of the combined observational samples Obs 1 and Obs2 
with the Millennium simulation samples Mill.AS and Mill.FS, using threshold density 6.0. 



old density 5.0. To avoid the inclusion to the combined obser- 
vational sample data of lower accuracy we used only super- 
clusters of the main sample with distances up to 520 /i"' Mpc. 
This combined sample Obsl contains 1151 superclusters, 592 
of them contain only 1 DF-cluster and were used in the calibra- 
tion of relative luminosities of superclusters (for this sample we 
getLo = 6.73e -I- 11). The other combined observational sample 
Obs2 was found using threshold density 6.0 and distance limit 
520 h^^ Mpc in all subsamples. 

The most striking feature of the Figure is the demonstra- 
tion of the presence of numerous very luminous superclusters 
in observational samples, and the absence of such systems in 
simulated samples. This difference between real and simulated 
supercluster richness is well seen using both richness criteria, 
the multiplicity and luminosity functions. For threshold den- 
sity 6.0 most luminous simulated superclusters have a relative 
luminosity of about 20 - 30 in terms of the mean luminosity of 
richness class 1 superclusters, most luminous superclusters of 
real samples have a relative luminosity about 100, i.e. they are 
about 3 times more luminous. The fraction of very luminous su- 



perclusters (relative luminosity 20 and above) is about ten times 
higher in real samples than in simulated samples. Similar differ- 
ences exist between the multiplicity functions of real and sim- 
ulated supercluster samples. For threshold density 6.0 the rich- 
est model superclusters have a multiplicity of about 20 whereas 
real superclusters have over 50. The number of Abell clusters in 
the richest superclusters is about 30; this difference in the num- 
ber of Abell and DF-clusters can be explained by differences in 
the number density of these cluster samples: the density of DF- 
clusters is about 2 times higher than that of Abell clusters. The 
differences between real and simulated samples are observed 
not only in the region of most luminous superclusters: over the 
whole richness scale the number of DF-clusters in simulated 
samples is smaller than in real superclusters. 

One more interesting observation: very luminous super- 
clusters are located in all subsamples (Northern and Southern 
regions of 2dFGRS, and in subregions of the SDSS DR4 sam- 
ple, if divided into 3 wedges of equal width). These subsamples 
have characterictic volumes of about 10 million cubic h^^ Mpc, 
whereas model samples of 10 times larger volume have no 
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extremely rich superclusters. Kolmogorov-Smirnov test shows 
that the probability that real and model distributions of super- 
cluster luminosities and richnesses are taken from the same par- 
ent distribution is 10"'" and 10"^, respectively. 

4. Luminous superclusters and inflation 

Superclusters of galaxies are formed by density perturbations 
of large scales. These perturbations evolve very slowly. As 
shown by Kofman & Shandarin ( 1988 1, the present structure 
on large scales is built-in already in the initial field of linear 
gravitational potential fluctuations. Actually they are remnants 
of the very early evolution and stem from the inflationary stage 
of the Universe (see Kofman et al. ll987> . The distribution of lu- 
minosities (or masses) of superclusters allows us to probe pro- 
cesses acting at these very early phases of the evolution of the 
Universe. 

One possible explanations for the large difference between 
the distribution of luminosities of real and simulated samples 
is the underestimate of the role of very large density perturba- 
tions. If this is the case then this means that our simulations 
have not yet reached a volume which can be treated as a fair 
sample of the Universe. In other words, a fair sample of the 
Universe has linear dimensions far in excess of 500 /i"' Mpc, 
used in simulations investigated in this Letter As shown by 
Power & Knebe J2006> . variations in the box size in a smaller 
box do not influence properties of Dark Matter haloes in cos- 
mological simulations. 

The other feasible explanation of the differences between 
models and reality may be the presence of some unknown pro- 
cesses in the very early Universe which give rise to the forma- 
tion of extremely luminous and massive superclusters. 

To date it is too early to make definite conclusions on the 
character of processes during inflation which may have caused 
the formation of very massive superclusters. Some aspects of 
this problem were recently studied by Saar et al. ("20061. They 
demonstrated that rich superclusters formed in places where 
large density waves combine in similar phases to generate high 
density peaks. The larger the wavelength of such phase syn- 
chronization, the higher the richness and mass of superclusters. 
The synchronization has properties of sound waves, where, in 
addition to the main frequency, overtones appear. In this con- 
text it is interesting to note that very rich superclusters have 
a tendency to form a quasi-regular network with characteristic 
scales 250 and 120 h^^ Mpc, as demonstrated by Broadhurst et 
al. (1990) and Einasto et al. ( 1994 1997b, 199^. 

The explanation of the physical origin of very massive su- 
perclusters is a challenge for theory. To get a more complete ob- 
servational picture of the phenomenon, large contiguous deep 
redshift surveys are needed. Only contiguous surveys allow 
to detect very massive superclusters. This is one reason why 
the continuation of the SDSS survey is so important, until the 
whole Northern hemisphere is covered by redshift data, as orig- 
inally planned. 
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